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Abstract 

This contribution describes the reorientation of the magnetization 
of thin ferromagnetic Heisenberg films as function of the temperature 
and/or an external field. Working in a rotating frame allows an ex- 
act treatment of the single-ion anisotropy when going to higher order 
Green's functions. Terms due to the exchange interaction are treated 
by a generalized Tyablicov (RPA) decoupling. 



1. Introduction 

New experimental results concerning the properties of thin magnetic films 
have provided an incentive for theoretical investigations. This contribu- 
tion deals with the reorientation of the magnetization of thin ferromagnetic 
Heisenberg films as function of the temperature or an external magnetic 
field. It is of particular importance to take into account collective exci- 
tations (magnons), which influence the magnetic properties of films more 
strongly than those of bulk magnets. Many-body Green's function theory 
(GET) is an appropriate tool to achieve this. In previous work, we investi- 
gated the reorientation transition by applying GET jHl^l, where a Tyablikov 
(RPA)-decoupling of the exchange interaction terms and a Anderson-Callen 
decoupling "3] of the single-ion anisotropy terms in the equation of motion for 
the lowest-order GE was used. By comparing with Quantum Monte Carlo 
(QMC) results |3] it was shown that this leads to rather good results for 
small anisotropies {K2 < O.IJ) when the magnetic field is in the direction 
of the anisotropy but the approximation is not as good when the field is ap- 
plied perpendicular to the anisotropy. A considerable simplification and an 
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improvement of the results concerning the reorientation is reported in Ref. 

where the Anderson- Callen decouphng is made in a frame which is ro- 
tated with respect to the original one and in which the magnetization is in 
the direction of the new z-axis. The reorientation angle is determined from 
the condition that the magnetization commutes with the Hamiltonian in the 
rotated frame. In this connection see also Ref. [U], who also apply the approx- 
imate Anderson-Callen decoupling in a rotated frame. We realized in Ref. 
[7j, generalizing findings of Ref.[H^, that when introducing higher-order GF's 
a decoupling of the anisotropy terms is not necessary. By taking advantage 
of relations between products of spin operators P, this leads to an automatic 
closure of the hierarchy of equations of motion with respect to the anisotropy 
terms. This procedure gives improved results (as compared to the approx- 
imate Anderson-Callen treatment) when applying a field in the direction of 
the anisotropy, which now can be large |7j, but there are difficulties (there 
were zeroes in the equation of motion matrix requiring a special treatment 
fOlE]; and numerical instabilities occurred) when trying to calculate the 
reorientation induced by a field perpendicular to the anisotropy. The new 
result of the present paper is that these difficulties can be overcome when 
working in a rotated frame, in analogy to Ref. [5]. In this way, we are able to 
describe the field-induced spin reorientation transition for spin S > 1 with 
an exact treatment of the single-ion anisotropy. The exchange interaction 
terms are still decoupled with generalized RPA. 

2. The Green's function formalism 

We investigate a spin Hamiltonian consisting of an isotropic Heisenberg ex- 
change interaction between nearest neighbour lattice sites with strength Jki, 
a second-order single-ion lattice anisotropy with strength i^2,fe vertical to the 
film {x — ?/)-plane and an external magnetic field B = (Sq , 0, B^): 

n = -\ E Jki(s,sr + sist) - Y.K2,k(sif - Y.(b^oSI + B^.si). (i) 

<kl> k k 

In this paper, we restrict ourselves to this Hamiltonian; we dealt with the 
dipole-dipole interaction in Ref. [2| and with exchange anisotropics in Ref. ^2] • 
Owing to the external i?Q-field, the magnetization vector, initially in the 
z-direction, will rotate by an angle 6 in the x — 2;-plane, pointing now in 
the 2;'-direction of a new frame {x',y',z'). As in Ref. [Sj, we shall do the 
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calculations in the primed system, in which the magnetization vector has the 
components (0,0, {S^ ')). The transformation between the frames is 

/ (5^) \ / cos^ sine \ / (5^ ') \ 

{sy) \= 1 (5^ ') • (2) 

V {S') J V -sine cose / V {S") J 

Because {S^ ') = {S^ ') = in the rotated frame, one needs only to calculate 
{S^ ') in order to find the components of the magnetization in the original 
frame provided the angle 6 is known. 

After transforming the Hamiltonian (Q) to the primed system the follow- 
ing Green's functions are needed 

^tj' = ii^t 'j ')) ) 
GSf+'- = {{{Srr-\2St'-l)St';S-')) . (3) 

The single-ion anisotropy is active for spins S > 1, and one needs the 
first Green's function plus those for n = 1,2,3, ... to treat films with S = 
1,3/2,2,.... 

In establishing the equations of motion, the exchange interaction terms 
are treated by a generalized Tyablikov (RPA)-decoupling, in which we do 
not break products of spin operators with equal indices 

{{{S^ rst '-,3; ')) ^ {{S! r){{st ')) + i^t r;Si ')) . (4) 

Note that in the rotated system, {S^ ') = 0; i.e. the second term vanishes. 
After applying the decoupling procedure and performing a Fourier transform 
to momentum space one obtains the following set of equations of motion. 

ujG+^- = 2{S' ') + {S' ')J{q - 71^)^+'- 

+ {B^Q sin e + B'^ cose)G+'- + i^2(l - | sin^ e)G^+'-, 
ujG'+^- = (6((5^ ')') - 25(S + 1)) - i J7k(6((^^ ')') - 2^(^ + 1))G+'- 

+Jq{S' + (fig^-sine + Bo^cose)G'^+'- 

+7^2(1 - 1 sin^ e) (2G'(^)'+'- - G'+'-) , 

+^7k(|5(S + 1) + {2S{S + 1) - 1)(S^ ') - ')') - 4((5^ 0')^+'- 

+Jq{S' ')G'(")'+'- + (5o^ sin e + cos e)G'(^)'+'- 



3 



fsin^ 



G(z? 



10((5" ') ) + 8((5^ ') ) - (65(5 + 1) - 5)((5^ ') ) 



-{AS{S + 1) - 1){S' ') - 5(5 + 1) 

+J7k(|5(5 + 1) + (2S(5 + 1) - \){S' ') + (35(5 + 1) - |)((5^ ')') 
-4((5" ')') -5((5^ 

+ Jg(5" + (5o" sin 9 + cos ^)G'(")'+'- 



I sin^ I 



(5) 



For a square lattice with lattice contant a = 1, the quantity 7k = 2(cosA;a; + 
cos ky), and q = 4:, the number of nearest neighbours. 

As Ref.jni we neglect all GF's not containing an equal number of S~ ' 
and 5+ ' operators. 

One observes that in the equations (0) the anisotropy terms do not yet 
lead to a closed system. This is, however, achieved by using formulas derived 
in Ref.p, which reduce products of spin operators by one order (!). One 
obtains 



for S = 1 : 




^ = 1(6-^+'- 


- + , 


for S = 3/2: 


Qizf + ,- 






for S = 2 : 












- gG 





(6) 

When inserting these relations into the system of equations (0) one sees that 
the resulting system of equations is closed. 

The equations of motion can be written in compact matrix notation 



(uj1-T)G = A. 



(7) 



The quantities F, G, and A can be read off from equation (0), where the 
non-symmetric matrix F is a (2 x 2), (3x3), (4x4) -matrix for spins 
5=1, 3/2, 2, respectively. The desired correlation vector corresponding to 
the Green's functions (jS)), 



C 



(5- '5+ ') 
(5- '^(^')"''(2^^ 1)^+ ') 



(8) 
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is obtained via the spectral theorem. Using the eigenvector method of Ref . [2] , 
one has for the components of the correlation vector C after a Fourier trans- 
form to configuration space 

dkCi(k) = — dhj. j dky ^ RijejkLkiAi (9) 

^ j,k,l=l 

(^ = 1,2,..,2^), 

where the integration is over the first Brillouin zone, and R(L) are matrices 
whose columns (rows) consist of the right (left) eigenvectors of the matrix F, 
and ejk = Sjk/ (e^"^^ — 1) is a diagonal matrix, in which Uj are the eigenvalues 
(j = 1, ..,25') of the F-matrix. 

It remains to derive an equation which determines the rotation angle. 
This is again done by using the approximation that the commutator of the 
magnetization with the Hamiltonian vanishes in the rotated system. This 
implies that the following Green's function is zero 

{{[Sr,H'];S-')) = 0. (10) 

Evaluating the commutator, one obtains 

{B^ cos e-B^ sin 6)0+^- - K2 sin 6 cos ^G-/'" = 0. (11) 

After applying the spectral theorem to this equation one obtains for the 
diagonal correlations the relation which determines the reorientation angle 
angle: 

{B^ cos 9 - B^ sin 9)C-^ - K2 sine cos 9C-'^ = . (12) 

This is the generalization of the angle condition of Refs. when treating 

the single-ion anisotropy exactly instead of applying the Anderson-Callen 
decoupling. Equation (fT^ together with the set of integral equations Q 
have to be solved self-consistently in order to obtain the magnetization {S^ ') 
and its moments in the rotated system together with the reorientation angle 
9. Applying the relations (0) yields the components of the magnetizations in 
the original system in which the magnetic reorientation is measured. 

3. Results 

We solve the equations of the last section with the curve following method de- 
scribed in detail in an appendix of Ref. 0. In Ref. [S], where the Anderson- 
Callen decoupling was used in the rotated frame, only cases with weak single- 
ion anisotropy were considered. An example was the case of spin 5 = 2 and 
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K2 = O.OIJ, which is a rather small anisotropy as appearing for 3d transi- 
tion metals. A very good agreement with QMC calculations of Ref. |1| was 
obtained. Because the anisotropy is small the result with the present theory 
is practically indistinguishable, whereas the results of Ref. [I] are rather bad 
for the reorientation, because the decoupling was not done in the rotated 
frame. For anisotropies in the rare earth region, which can be of the order 



Compare: K2(exact) o , A.C.-decoupling ^ , QMC • 
S=2, J=1 00, K2=20 T=1 00 
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Figure 1: Normalized magnetizations (S^) / S and {S^) / S and the reorien- 
tation angle G for a spin S=2 Heisenberg monolayer as function of the ex- 
ternal field are shown: QMC ^ (solid circles), Anderson- Callen decoupling 
[S] (triangles) , present theory (open circles). 

of the exchange interaction, the approximate theory of Ref. j3] should break 
down, and one expects the present theory to be superior. Surprisingly, the 
Anderson-Callen decoupling in the rotated frame still yields excellent results, 
when compared with the present theory and QMC results from Ref. ^ for 
anisotropies up to K2 < 0.2 J. This can be seen in Fig. 1 for the magnetic re- 
orientation induced by the transversal i?^-field for the case of K2 = 0.2 J and 
T = 100. The result of both Green's function theories are nearly identical 
to each other and deviate only slightly from the exact ( within the statistical 
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error) Quantum Monte Carlo results. 
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Figure 2: Normalized magnetizations {S^)/S and {S^)/S and the reorien- 
tation angle G/(7r/2) for a spin S=2 Heisenberg monolayer as function of a 
transversal field are shown: Anderson-Callen decoupling (dotted lines) 
and the present theory (solid lines) for K2 = 0.5J (T/J = 4.9). 



The reason that the magnetization curves as function of the external 
field of both theories are very close for T=100 is that with the present 
choice of the parameters the magnetizations at this temperature are still 
very close to each other. Both theories deviate from each other with increas- 
ing temperature, until close to the reorientation temperature the deviation 
is maximal. Therefore one should observe deviations between the theories 
in this temperature range. The differences should also increase with increas- 
ing anisotropy strength K2. Therefore we have calculated the field-induced 
reorientation for a temperature somewhat below the reorientation temper- 
ature and for a large anisotropy K2 = 0.5 J (T/J = 4.9), for a Heisenberg 
monolayer with spin = 2. The result is shown in Fig. 2. In this case the 
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Anderson-Callen (A.C.) decoupling along the lines of [3| leads to a discon- 
tinuous transition from a certain angle {6/ {71/2) ^ 0.6) to full reorientation 
(6'/(7r/2) = 1). whereas the reorientation transition is continuous when the 
single-ion anisotropy is treated exactly. Such discrete transitions are also re- 
ported in Ref. [HI in a treatment which is very similar to that of Ref. [5j. The 
reason why this is not observed in the last reference is that there only very 
small anisotropies were considered. We attribute the discontinuous transi- 
tion to the approximate Anderson-Callen decoupling, which is not justified 
when going to large anisotropies. The difference between the corresponding 
reorientation fields, B^, increases with increasing anisotropy. For the present 
case it is: S^''^'""* - B^-^- ~ 11 (for K2 = 0.5 J). 

Unfortunately, we cannot say how accurate the present model is because 
we have no QMC calculations available for large anisotropies. The most un- 
controlled approximation consists in the generalised RPA decoupling of the 
exchange interaction terms of the higher-order GF, eqn.(j3]) . Previous cal- 
culations PHI have shown (by comparing with QMC) that RPA is a good 
approximation for a Heisenberg model (no anisotropy) with a field perpen- 
dicular to the film plane. For improving the the present approach for a field 
in transversal direction one may try to extend the present theory by apply- 
ing the procedure of jT3] which goes beyond the RPA with respect to the 
exchange interaction terms. It is possible by applying formulas from Ref.jH] 
to treat also the fourth-order anisotropy term, — J2i -^4,j('S'f )^, exactly. Work 
in this direction is in progress, as well as the generalization to multilayers. 
We mention that the reorientation of the magnetization of a ferromagnetic 
film induced by the coupling to an antiferromagnetic film can also be treated 
with the GF formalism |15j . 
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